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ABSTRACT 

Climatic dynamic topography variations in the Alaska gyre during the period 1968-1990 are described with 
an objective analysis of more than 12 000 STD and XBT stations, and COADS wind stress data. Interannual 
dynamic height and SST variations were correlated and were consistent with recently described large-scale 
climatic shifts in the North Pacific. The gyre was centered more to the east, circulation appeared stronger, and 
SST was lower during the early to mid-1970s than during the 1980s. The Aleutian low (NP and PNA indices) 
intensified during the interim, but the response did not appear as a gyre spinup. Instead, the associated wind 
stress anomalies forced a slowly varying dynamic height anomaly across the eastern and northern part of the 
gyre through Ekman convergence, which had the effect of displacing the gyre’s low somewhat to the WSW in 
the 1980s. The wind curl spectrum was white, and the slow oceanic response was modeled as stochastic-forced 
climate variability with a simple first-order Markov autoregression process. Forcing was assumed to be Ekman 
pumping of the pycnocline, and the damping coefficient was estimated from the data to be ~ 1 yr . A hindcast 
with observed winds gave estimated dynamic height patterns similar to those observed, with a canonical cor- 
relation of 0.79 at 99% confidence. This response was weak in the western half of the gyre, where slow baroclinic 
variability may have been influenced by long Rossby wave propagation. A simple autoregression simulation 
using artificial white noise forcing shows the evolution of decadal variations similar in nature to those observed. 
This result, along with the low frequency correlation between dynamic height and SST, suggests that the upper- 
ocean climatic variability in this region is primarily wind forced. 


1. Introduction 

Oceanic variability on decadal (sub-ENSO) time- 
scales can influence weather, fishery resources, other 
socioeconomic factors, and obscure climate trends, 
which many believe are the consequence of anthropo- 
genic influence on the environment. Wunsch (1992) 
argued that interdecadal oceanic variability exists as a 
natural outcome of the interactions between turbulent 
fluids such as the ocean and atmosphere, but noted that 
the present oceanographic database is generally inad- 
equate for a reliable description or understanding. 
There has been progress, such as in the comparatively 
data-rich North Atlantic, where Levitus (1989) docu- 
mented large-scale middepth circulation changes be- 
tween decades. However, this was based on the differ- 
ence between pentadal averages and did not offer clues 
on the variability timescales such as could be obtained 
with a continuous time series. Sea surface temperature 
(SST) is the most extensively observed oceanic vari- 
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able for climate analysis. Although encumbered with 
calibration problems earlier in the record, these data 
have been applied by numerous authors in correlation 
studies to assess the effect of SST anomalies on the 
atmosphere and vice versa (Kushnir 1994, to cite a re- 
cent example). 

SST is strongly influenced by air- sea interaction 
and thus is not the best indicator of ocean circulation 
variations. Dynamic topography is more appropriate. 
Spacebome altimetry is expected to accumulate long 
continuous global records suitable for climate studies 
in the future (Koblinsky et al. 1992). In the mean- 
time, it is possible to generate dynamic topography 
analyses in regions relatively densely sampled by 
STD and XBT observations. This study applies such 
an analysis to the Alaska gyre, where interdecadal 
SST variations have been documented and a limited 
number of conventional oceanographic surveys have 
shown significant transport variations in the major 
currents ( see below) . A large XBT data resource ex- 
ists, and an analysis is presented here of the available 
hydrographic stations spanning the years 1968 
through mid- 1990; 77% XBTs and 23% STD/Nan- 
sen stations (see the appendix). From these, a 22.5 
year long horizontally gridded field of 0-450 decibar 
(dbar) dynamic height has been computed. (1 dbar 
= 10 4 Pa). The focus here is to examine gyre-scale 
dynamic topography variations during recent dec- 
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ades, the relation to climate indices, and the dynam- 
ics of wind forcing. 

Of particular interest is the low-ffequency oceanic 
response to high-frequency wind forcing. The problem 
lends itself to application of stochastic forcing models 
(Hasselmann 1977), which can be evaluated with ex- 
tensive dynamic height and wind measurements. The 
latter consist of monthly wind stress fields from Cou- 
pled Ocean- Atmosphere Data Set (COADS). It is ev- 
ident in these data and other studies of model winds 
(Chave et al. 1991; Large et ai. 1991) that the wind 
curl spectrum is white and the forcing appears stochas- 
tic (section 5). Hasselmann ( 1977) introduced the no- 
tion that slow climatic variability can be explained as 
the integral response to continuous random short-term 
impulses. A large frequency separation between forc- 
ing and response is an essential element of this theory. 
Statistically stationary climate variability is maintained 
with linear damping feedback and the problem ex- 
pressed as a simple first-order Markov process. The 
theory has been applied to explain SST variability in 
response to random air- sea fluxes (Frankignoul and 
Hasselmann 1977; Frankignoul 1981; Frankignoul and 
Reynolds 1983). The Markov model is explored here 
to explain large-scale dynamic height variations as the 
response to stochastic vertical pycnocline displacement 
driven by Ekman pumping. Reed (1984) suggested 
time-integrated Ekman pumping might explain ob- 
served transport variations in the eastern Alaska Gyre 
between 1980 and 1981. Other studies of stochastic 
wind forcing have focused on the mesoscale response, 
including remote forcing, using more sophisticated 
models (Frankignoul and Muller 1979; Muller and 
Frankignoul 1981; Treguier and Hua 1987; Brink 1989; 
Samelson 1989, 1990; Cummins 1991; Large et al. 
1991; Cummins and Freeland 1993). Here, the data 
duration and spatial extent allow the climatic response 
at lower frequencies and wavenumbers to be explored. 
Frankignoul and Muller (1979) showed that long 
Rossby wave resonance is possible at periods longer 
than the minimum Rossby wave period (T min — 1.9 
years at 55° latitude and internal radius L D ~ 16 km). 
Thus, we can anticipate a slow oceanic response to the 
low frequency and wavenumber modes in the atmo- 
sphere. 

For the present application, the Markov model is ex- 
pressed as 

dn 

— = -aW £ - \rj, ( 1 ) 

where 77 is the dynamic height, W E the stochastic Ek- 
man pumping velocity, a a scaling coefficient, and \ 
the damping or feedback coefficient. Here \ > 0 in- 
duces a negative growth tendency to a positive anomaly 
and vice versa; \ has units of inverse time, and it is 
often convenient to express it in terms of a damping 
timescale At this point we will defer discussion 


on what physical processes play this role, as for ex- 
ample the atmospheric feedback proposed for stochas- 
tically forced SST anomalies by Frankignoul and Has- 
selmann (1977). The important features of the 
two-time-scale climate system can be addressed inde- 
pendently of the details of the upper-ocean model 
(Frankignoul and Hasselmann 1977), and it is only 
necessary to accept that some negative feedback mech- 
anism exists. To test this model, the observed fields 77 
and W E were used to empirically estimate the coeffi- 
cients a and \ at each grid point. These exhibit spatial 
patterns indicating that the wind-forced response was 
not uniform. The Markov model is then used with these 
coefficients to hindcast dynamic height from the ob- 
served winds. This hindcast field shows very good 
agreement with the observed field in both temporal and 
spatial pattern (section 5). 

2 . Regional ocean circulation and climate variability 

The Alaska gyre is a subpolar cyclonic circulation 
in the northeast Pacific (Fig. 1) associated with the 
Aleutian low, the large-scale winter low pressure sys- 
tem centered over the Aleutian Islands, some 2000 km 
westward. The primary currents consist of a broad east- 
ern boundary current flowing north, condensing into a 
narrow western boundary flow in the apex of the gyre 
and proceeding WSW along the Aleutian Peninsula as 
the Alaskan Stream. Geostrophic transports relative to 
1500 dbar in the Alaskan Stream are nominally about 
12-18 Sv (Sv » 10 6 m 3 s~') (Musgrave et al. 1992). 
Transports as small as 2 Sv and as large as 20 Sv have 
been observed from intermittent surveys (Reed et al., 
1980; Royer 1981). Reed et al. (1980) and Musgrave 



Fig. 1 . Mean 0-450 dbar dynamic topography in the Alaska gyre 
as computed in this analysis (the appendix). Inset: the study area 
location (shaded) in the northeast Pacific. The letter “A” denotes the 
approximate mean center of the wintertime Aleutian low of surface 
atmospheric pressure. 
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et al. (1992) found that this variability was uncorre- 
lated with seasonal Sverdrup transport estimated from 
integrated wind stress curl. Reed (1984) speculated 
that vertical pycnocline displacements, associated with 
shifting patterns of wind stress curl, induced east- west 
lateral displacements of the northward transport feed- 
ing the gyre. The importance of local wind stress curl 
forcing to episodic ventilation of the upper pycnocline 
was suggested by Van Scoy et al. (1991). Apparent 
lateral displacements were also evident in the early 
charts of Dodimead ( 1963). Royer and Emery ( 1987) 
suggested that these could be associated with meridi- 
onal shifts in the North Pacific drift interacting baro- 
tropically with bottom topography. Musgrave et al. 

( 1992) argue that the apparent transport variations may 
be an artifact of undersampling the field in the presence 
of large eddies. The Geosat altimeter provided high res- 
olution sampling from 1986 to 1989. Kelly et al. 
(1993) described a principal mode with a large-scale 
sea level increase across the NE part of the gyre. Bhas- 
karan et al. (1993) found modes of interannual varia- 
tion that correlate to climate indices such as the North 
Pacific pressure index (NP) and Southern Oscillation 
index (SOI). 

SST variations of 1-3 decades duration have been 
identified in the gyre by Royer (1989). Correlation 
with Alaskan coastal air temperature records allowed 
extrapolation back to ~1900, wherein the period ap- 
pears to have modulated from as short as 9 to as long 
as 23 years, and a nominal value of 16 years. The quasi- 
cyclic behavior has led to speculation of a link to har- 
monic astronomical forcing with the 18.6 year nodal 
tide, wherein weak meridional barotropic currents in 
the presence of the mean SST gradient could induce 
SST fluctuations on the order of 2°C (Royer 1993). In 
light of this, it is instructive to examine alternative pro- 
cesses to determine whether these SST variations are 
related to large-scale variations in the gyre circulation 
and the atmospheric forcing rather than very low fre- 
quency tides. Figure 2a shows the SST time series (2- 
year running mean) from the 5° latitude /longitude box 
in the east-central part of the gyre (centered at 55°N, 
140°W, after Royer 1989). Superimposed is the first 
empirical orthogonal function (EOF) of the gyre dy- 
namic height from the present study ( see section 4 for 
details). It is immediately apparent that these indices 
share similar variations on interannual to interdecadal 
timescales (correlation r = .74), suggesting that SST 
is closely linked to the dynamic height field and related 
processes governing the gyre circulation. 

Trenberth and Hurrell ( 1994) and Miller et al. 
( 1994) describe the climatic shift in the North Pacific, 
which was manifested rather abruptly in the mid-1970s 
and persisted for about the next decade. The symptoms 
included deepening of the Aleutian low sea level pres- 
sure (SLP), southward displacement of Pacific storm 
tracks, lower SST, and deeper mixed layer across the 
central North Pacific and higher SST along the North 



Fig. 2. (a) Light line: Time series of monthly SST anomaly in the 
5° X 5° square centered at 55°N, 140°W, with a two-year running 
mean filter applied (monthly 5° X 5° SST data courtesy of D. Cayan. 
Scripps Institution of Oceanography). Dark line: Superimposed time 
series of the first dynamic height EOF from the present analysis. The 
similarity suggests that both indices reflect similar oceanic response 
on decadal timescales. The horizontal dashed line is the offset for the 
SST mean during the time interval from 1968 to 1990 and serves as 
the abscissa for the dynamic height EOF curve, (b) Solid line: The 
NP index (Trenberth and Hurrel 1994) shown negative. Dashed line: 
The PNA index (Wallace and Gutzler 1981). Both are 5-year running 
mean monthly anomalies and indicate minima in the early to mid- 
1970s and increase during the 1980s. An increase as plotted signifies 
deepening of the Aleutian low. [The signal is more robust in the NP, 
which is an area-average surface pressure index than the PNA, which 
is an upper air pressure teleconnection index. (Data courtesy of K. 
Trenberth, NCAR and M. Wallace, University of Washington).] 


American coast and Alaska, along with changes in a 
host of other indices. The Pacific North America 
(PNA) mode (Wallace and Gutzler 1981) and the 
North Pacific (NP) index (Trenberth and Hurrell 1994) 
show trends indicating a deepening Aleutian low over 
the decades of the 1970s and 1980s (Fig. 2b). Tren- 
berth ( 1990) showed that the 1970s climate transition 
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was associated with warm surface temperature anom- 
alies over Alaska and cooler SST in the central North 
Pacific. Wallace et al. ( 1990) noted that the first EOF 
North Pacific SST pattern is robust, having been evi- 
dent in several previous studies using various analysis 
techniques, domains, and periods. The pattern reflects 
an elongated region dominating the central North Pa- 
cific, with a reversal of sign in the northeastern basin 
and the Gulf of Alaska, thus indicating an opposite 
phase in the low-frequency interannual SST variability 
between the Alaska gyre and the central North Pacific. 
Their temporal expansion of that mode pointed to a 
cooler Alaska gyre SST in the 1970s and a wanner one 
in the 1980s, as indicated here in Fig. 2a. They also 
noted that winter SST tendency (SST change from 
early to late winter) correlations with 500-mb height 
anomalies favor large-scale atmospheric forcing of the 
ocean. In the Gulf of Alaska, these correlations indi- 
cated a winter SST warming tendency associated with 
a stronger Aleutian low (+PNA). 

3. Data processing and mean fields 

The set of 9382 XBT and 2769 STD stations ob- 
tained from the National Oceanographic Data Center 
( see the appendix ) were used to compute 0-450 dbar 
dynamic height, t/(x, y, t), within the boundaries 50°N 
latitude, 160°W longitude and the coast. Converting 
XBT temperature profiles to dynamic height is de- 
scribed by Lagerloef (1994). An rms error of <0.03 
dynamic meters (dyn m) was obtained using a regres- 
sion model relating empirical dynamic height modes to 
integrated temperature profiles. This empirical halo- 
steric adjustment was a significant improvement over 
the conventional temperature salinity correlation 
method in the Alaska gyre. Even with such adjustment, 
these dynamic heights reflect the integrated heat con- 
tent of the upper ocean to a large extent, and the linkage 
with the SST anomalies is anticipated (Fig. 2a). Dy- 
namic heights were then interpolated to a horizontal ( 1° 
lat X 1° long) and temporal grid by an objective anal- 
ysis using decorrelation scales of 200 km and 90 days 
(the appendix). Figure 1 shows the mean dynamic 
height relative to 450 dbar reflecting the general cy- 
clonic circulation of the gyre. 

Wind stress data were derived from the COADS 
monthly mean 2 deg x 2 deg pseudostress (WU and 
WV) product available from the National Center for 
Atmospheric Research. These were converted to wind 
stress components r t and r v by reducing WU and WV 
to equivalent wind speed terms U, V and W = (U 2 
+ V 2 ) 1/2 and applying the Large and Pond ( 1981 ) for- 
mulation. The stress components were then interpo- 
lated to the time and space grid used for the ocean 
dynamic height data, and applying additional smooth- 
ing of 90-day and 1000-km decorrelation scales. From 
these, the mean stress, Ekman transport components 
and Ekman pumping velocity were obtained from stan- 
dard formulas (Fig. 3). 


The Ekman transport (Fig. 3b) reveals a highly di- 
vergent field with an origin of flux in the northern part 
of the gyre. Flux is ~ 1 Sv southward out of the gyre 
across the southern boundary at 50°N, associated with 
the prevailing westerlies. The Ekman pumping veloc- 
ities W E (Fig. 3c) are maximum in the northern regions 
where divergence is concentrated and are comparable 
to or exceeding published values for equivalent wind 
stress curl, pfW E , in the North Pacific (Large et al. 
1991; Chave et al. 1991; Rienecker and Ehert 1988). 
The area integral of W E yields a net vertical transport 
comparable to the net Ekman flux out of the gyre. There 
is an evident process of localized mean Ekman pump- 
ing in the north extracting mass from the interior geo- 
strophic flow and conveying it southward in the surface 
Ekman flux. 

4. Low-frequency variations 

Taking the conventional approach to studying inter- 
annual variability, data presented henceforth are anom- 
alies from seasonally averaged values. As noted above 
in Fig. 2a, 77 EOF1 is an expression of similar inter- 
decadal variability previously identified in SST. This 
mode (shown on expanded scale in Fig. 4a) is domi- 
nated by an ~20-yr period punctuated by shorter mul- 
tiyear variations. It accounts for 29% of the dynamic 
height variance, which is concentrated in the northern 
and eastern portions of the gyre. The 77 EOF2 ( 8% vari- 
ance) reflects a large feature in the north and northwest 
part of the gyre, the negative phase culminating in the 
mid- to late 1970s and the positive persisting after 
mode 1 has declined in the late 1980s (Fig. 4). It 
should be noted that these fields are considerably 
smoothed by the filtering effect of EOFs. There is con- 
siderable mesoscale variability in the higher modes of 
these seasonal anomalies, as well as a seasonal cycle, 
which are not being addressed here. Modes 1 and 2 
encompass the primary low-frequency decadal varia- 
tions that are of interest. 

A chronology of these variations is clarified with an 
overlapping sequence of 4-year mean dynamic topog- 
raphy anomalies in Fig. 5. Throughout the early and 
mid-1970s a large-scale negative anomaly persisted 
across the basin with a slow migration from east to 
north to west. Toward the late 1970s and into the early 
1980s a positive anomaly accumulated in the southeast 
and spread north and westward. By the mid-1980s the 
positive anomaly dominated the basin much the way 
the negative anomaly had done a decade before, and 
toward the late 1980s negative anomalies appeared in 
the south and east. Thus, the 1970s encompass a cool 
episode, with low dynamic height and SST, while the 
1980s were warm, with large-scale positive dynamic 
height anomalies. Minimum and maximum conditions 
are represented by the 4-year means marked by 1976 
and 1984, respectively. The sum of these anomalies 
with the mean dynamic topography illustrates the dif- 
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Fig. 3. (a) Mean annual surface wind stress adapted from 
monthly COADS pseudostress and interpolated to the present 
grid. The pattern reflects the mean annual atmospheric cir- 
culation associated with the Aleutian low. (b) Mean annual 
surface Ekman transport field computed from the mean stress 
field. Note the acute divergence concentrated in the northern 
region and the net surface transport of about 1 Sv southward 
across 50°N. (c) Mean annual vertical Ekman pumping veloc- 
ity computed from the divergence of the Ekman transport 
field. 


ferent flow conditions (Fig. 6). In the 4-year mean 
marked by 1976 the gyre was relatively intense with a 
well-defined center and strong Alaskan Stream. The pe- 
riod marked by 1984 showed a relatively relaxed cir- 
culation and the center displaced southwestward. These 
changes took place as the Aleutian low pressure system 
was intensifying between the 1970s and 1980s (Fig. 
2b). The ocean data indicate that the gyre filled rather 
than deepened during this period, counter to the con- 
ventional view that the strengthening Aleutian low 
would spin up the gyre. Warmer conditions in the 
northeast Gulf of Alaska in conjunction with a stronger 
Aleutian low are also consistent with patterns seen by 
Wallace et al. ( 1990) and others. 

The evolving wind stress anomaly field is depicted 
with a sequence of overlapping 4-year means in Fig. 7 
for comparison with Fig. 5. The early to mid-1970s 
were dominated by north, northwest, and westerly wind 
anomalies. In the latter decade the wind anomalies 
shifted to east and southeasterlies as the Aleutian low 
intensified. Ekman pumping was evidently the primary 


wind forcing mechanism according to these patterns. 
Generally a southward and divergent Ekman transport 
anomaly persisted in the 1970s, with positive W E in- 
ducing upward pycnocline displacement and reduced 
dynamic height. The converse was true in the 1980s: 
generally north and east Ekman transport induced con- 
vergence and accumulation of warmer buoyant waters 
in the upper layer. Upweiling anomalies may not en- 
tirely account for the SST variations between the dec- 
ades because a shallow isothermal or inversion layer 
often exists in the upper halocline of the northern Gulf 
of Alaska. SST may have been influenced by the sur- 
face Ekman transport anomalies, which advected 
cooler water from the north in the 1970s and wanner 
water from the south in the 1980s. A canonical corre- 
lation analysis (CCA) was applied here to determine 
modes maximizing the correlation between 77 and low- 
pass filtered (two-year running mean) r, using three 
EOFs input from each. [See Barnett and Preisendorfer 
( 1987) or Graham et al. ( 1987) for a description of the 
analysis.] The r CCA1 spatial amplitude depicted a 
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Fig. 4 . EOFs of dynamic height (77) anomalies from seasonal 
means, (a) Time series of 77 EOF1 (dark) and 77 EOF2 (light), con- 
taining 29 % and 8% of the variance, respectively, (b) 77 EOF! spatial 
pattern showing the positive anomaly across the east and north 
regions of the gyre, (c) 77 EOF2 spatial pattern showing a major fea- 
ture in the north and northwest part of the gyre (dashed contours are 
negative). 


NW-SE oriented wind pattern, with a zone of maxi- 
mum stress and minimal curl dominating the central 
portion of the gyre, generally clockwise curl to the east 
and counterclockwise to the west (Fig. 8 a). The r 
CCA 1 and 77 CCA1 time series amplitudes (Fig. 8 b) 
showed the same negative phase during the 1970s and 
positive in the 1980s (correlation r — .69). The cor- 
responding 77 CCA1 spatial amplitude (not shown) was 
nearly identical to 77 EOF1 (see Fig. 4b). 

5. Transient response to stochastic forcing 

An inspection of Fig. 8 b indicates a slight lag be- 
tween r CCA 1 and 77 CCA 1 time series, even though 
they represent signals with the highest in-phase corre- 
lation computed by the CCA. The delayed response 
was examined further with CCA computed at various 
lags. The 77 input was limited to EOF1 in this case be- 
cause it is the dominant signal to be explained and be- 
cause 77 EOF2 appears to represent a propagating fea- 
ture ( see section 6 ) that otherwise influences the lagged 
CCA in a manner not representing direct forcing. Max- 
imum canonical correlation of r « 0.8 was found with 
r leading 77 by — 1 year. Secondary peaks with 77 lead- 
ing r by ~3 years (r » 0.7) and r leading 77 by ~4 
years (r « 0.58) are related to an evident 3-4 year 
variation, which will be the subject of another inves- 
tigation and not discussed further here. The primary 
peak at 1 -year lag gives a strong indication that the 
dynamic height variations are a delayed response to 
wind forcing. 

To study this transient response with the Markov 
model, Eq. ( 1 ) , it was necessary to compute the Ekman 
pumping velocity W E . This was done using a finite- 
difference estimate of Ekman transport divergence. The 
wind stress data were unfiltered, except for the 90-day 
decorrelation scale applied with the objective analysis, 
and retained an unknown level of spurious measure- 
ment noise. To minimize the influence of this error on 
subsequent estimates of the curl, the stress components 
r j and r v were synthesized from their first four EOF 
modes (89% variance retained) prior to computing the 
Ekman terms. (The truncation criteria followed the 
usual practice of retaining modes whose Eigenvalues 
remained above those from EOFs of white noise.) The 
resultant W E anomalies are interpolated to the same grid 
as the other variables except that points within one 
grid space of the coast were lost due to the finite dif- 
ferencing. 

It is important to briefly discuss the predominant 
characteristics of the variability. Here W E EOF1 shows 
a concentric pattern centered over the gyre, while W E 
EOF2 shows broad downwelling across the northeast 
and upwelling to the west (Fig. 9). In the context of 
the Ekman pumping hypothesis, the W E EOF2 rather 
than the W E EOF l pattern is more consistent with the 
t CCA1 and 77 EOF 1 patterns. The two-year running 
mean of the W E EOF amplitude time series depicts the 
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Fig. 5. Sequence of four-year running means of tj anomalies centered at the start of the year 
indicated in each frame (e.g., 1970 represents the average of years 1968, 1968, 1970, 1971). Note 
the apparent counterclockwise progression of a large-scale negative anomaly around the gyre in 
the 1 970s. followed by the large-scale positive anomaly in the 1 980s (dashed contours are negative). 


low-frequency components of these modes (Fig. 10). 
W E EOF1 displays variations of ~5 yr periods, while 
W E EOF2 shows a low-frequency trend similar to the 
increasing PNA index from the 1970s to the 1980s (cf. 
Fig. 2b). Thus W E EOF2 appears to carry much of the 
subtle low-frequency variability associated with the 
longer-scale climate variations. In contrast, the unfiltered 


W E EOF time series are clearly dominated by much 
shorter timescales (Fig. 11), where low frequencies 
have much smaller amplitude than the shorter timescale 
atmospheric noise. Composite spectra of tj and W E (Fig. 
12) demonstrate the redness of oceanic variability com- 
pared to the white W E spectrum. This timescale separa- 
tion is characteristic of stochastic climate forcing. 
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Fig. 6. Four-year mean rj anomalies added to the mean annual 
dynamic height field for the periods marked 1976 and 1984 in Ftg. 
5. (a) In the four years marked by 1976 the gyre was relatively in- 
tensified with a well-defined center and strong Alaskan Stream, (b) 
The period marked by 1984 showed relatively relaxed circulation and 
the center displaced southwestward. 


Ekman convergence results in a — W Ej which should 
induce a positive drj/dt through pycnocline motion. 
This implies that the low-frequency variations should 
be out of phase by tt /2 and 77 should lag — W E , whereas 
the CCA analysis above showed a high degree of cor- 
relation nearly in-phase. This can be reconciled by ex- 
amining the stochastic forcing model ( 1 ) at each end 
of the spectrum. At high frequencies u> X, 

dri 

f t ~-«W E ( 2 ) 

and the response 77 is an accumulation of random im- 
pulses (random walk); 77 is theoretically in quadrature 
with - W E and the spectral slope is w ~ 2 . At frequencies 
w < X, 


X 77 » ~aW E , ( 3 ) 

and the response is approximately in phase with -W E . 
The low-frequency end of the spectrum becomes white 
and the variance statistically stationary (Hasselmann 
1977). The subrange of frequencies u> ~ \ marks the 
transition from red to white spectra and quadrature to 
in-phase response. It will be shown below that X 
** 1 yrj' in the gyre, with a corresponding period T 
» 27 rX 1 » 6 yr. Thus the subrange u ~ X encom- 
passes the interannual to interdecadal band, which 
dominates the observed variability. Note that Tis about 
three times larger than the minimum Rossby wave pe- 
T mm » 1.9 yr, which is the high-frequency dy- 
namic limit of the baroclinic response. The effective 
timescale of the low-pass filtering process by the ocean 
is defined by u> ~ X. Figure 12 illustrates an idealized 
frequency spectrum F(u) * (w 2 + X 2 ) -1 (see Fran- 
kignoul 1981), with X = 1 yr"', superimposed on 77 . 
As X decreases, the transition portion of the spectrum 
shifts to lower frequencies and the decorrelation time- 
scale of the ocean response increases. It can be shown 
that the theoretical phase angle between 77 and W E in 
Eq. ( 1 ) is 9 = tan ~ ! (w/X) , which tends toward zero as 
w 0. The phase lag in units of time for a given fre- 
quency is 9/u, which asymptotically approaches X -1 , 
or ~ 1 yr, as u -* 0; consistent with the ~ 1 yr lag CCA 
given above. 

In discrete form, ( 1 ) is expressed as 

rj(t -f- 1) - 77 (f) 

— = —aW E (t) - Xi 7 (r) (4) 

or 

7]{t + 1) = <f>ri(t) - <xW E (t)dt, (5) 

where $ = 1 — \dt, and (t + 1) and (r) indicate dis- 
crete time steps separated by dt. When W E is uncorre- 
lated white noise, this is a simple first-order autore- 
gression (AR1 ) process (Wei 1990); 0 is the AR1 co- 
efficient equivalent to the lag 1 autocorrelation 
coefficient. Ideally, the lag interval dt should be the 
natural decorrelation timescale for W E . This is difficult 
to properly determine here because the original 
COADS wind data are monthly means and were sub- 
sequently smoothed further in the objective analysis 
(the appendix), where the sample interval was arbi- 
trarily selected ( dt - 1/16 year) and not prescribed by 
the natural decorrelation timescale. Estimating <b 
through the autocorrelation function was therefore not 
satisfactory in this case, and a different approach was 
taken. 

By recursive substitution in (5), it can.be shown that 
77 at any time step (t + 1) is determined by the sum- 
mation over N previous time steps: 

77(r + 1) = <f> N+i -q{t - N) - ocdt £ rw E {t - n). 


(6) 
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Fig. 7. Sequence of four-year running mean wind stress (t) anomalies centered at the start of 
the year indicated in each frame, as in Fig. 5. The early to mid-1970s were dominated by north, 
northwest, and westerly wind anomalies. In the latter decade the wind anomalies shifted to east 
and southeasterlies as the Aleutian low intensified. 


Note that with 0 < <{> < 1, the first term on the right 
vanishes for large N and the function becomes inde- 
pendent of the initial state at r - N. In any case, the 
initial condition t]( t — N) — 0 can be assumed with no 
loss of generality and the term can be ignored. In the 
remaining summation terms, a is seen as a scaling co- 
efficient between W E and 77, which simply determines 
the amplitude of the response function. The weights 


serve as a filter. For <f> -*• 0 (\ -*• 1 Idt) the weights tail 
off sharply; the function becomes independent of its 
history and approaches a likeness of W E . Conversely, 
as<£-*l(\-+ 0)the function is highly dependent on 
its past and approaches an ordinary random walk with 
inverse square spectral slope. 

Coefficients 4> and a can be estimated from (5) by 
multiple linear regression in the form z - <f>x + ay. 
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Fig. 8. (a) Spatial pattern for r CCA 1 , the first canonical com 
tion analysis (CCA) mode between r (two-year running mean) ; 
rj. A zone of maximum stress and weak curl dominated the cen 
portion of the gyre, with generally clockwise curl to the east ; 
counterclockwise to the west (in the positive phase). This repress 
. w,n ^ t ^ s anomaly mode pattern having the greatest correlat 
with p. ft) The time series amplitudes of r CCA1 (dark) and n CC 
(light), showing the low anomaly in the 1970 s and high in the 191 
as seen in rj EOF 1. 


where 2 - r?(f + 1 ), * = 77(1) and y = - W E (t)dt . 
Subsequently applying these in (5) to hindcast 77 , it 
became apparent that a was seriously underestimated 
and slightly overestimated, owing to the high 77 au- 
tocorrelation and weak 77 versus W E cross correlation. 
An iterative procedure was devised to find the best fit 
between the observed 77 and the hindcast computed 
from W E , as follows. Step 1: An initial guess of 0 was 
taken from the above multiple regression. Step 2: A 
hindcast estimate Tj e was computed with 

le{t + 1) = <f>Tj e (t) - W E (t)dt , ( 7 ) 

which is equivalent to (5) except for omitting the scal- 
ing coefficient a. Step 3: a was then estimated by linear 


regression of tj - arj e , which yielded a larger value 
than the initial multiple regression. Step 4: A new es- 

^ ^ was obtained b Y bnear regression of 
{aWzi^di + rj(t + 1)} = This reduced d> 

slightly as a increased. Steps 2-4 were repeated until 
the rms difference between 77 and art, ceased to change 
significantly. Resulting values of (f> ranged from 0 84 
to 0.98 with a median of 0.94. The corresponding 
damping coefficient time scale = (1 - <f>)/d t fell 
in the range 0.5 -2.5 yr with a median value of ~ 1 y r . 

The scaling coefficient a was on the order of 10 “ 2 - 
10 , or ~g'/g (where g' is reduced gravity) as 
would be expected from dynamic scaling arguments. 
However, the hindcast amplitudes arj e tended to be 
about a factor of 2 smaller than the observed 77 time 
series. From (3 ), a represents the low-frequency re- 
gression between — W E and \tj. The spatial pattern was 
found to be nonuniform (Fig. 13), with the region of 
strongest response concentrated in the eastern half of 
the gyre. This eastward bias is consistent with 77 EOF1 
and the implicit Ekman pumping in W E EOF2. A large 
region enclosed by a » 0 signifies that the local wind 
forced response is weak through much of the western 
S5* sy re - This is the region dominated by W E 
bUrl, indicating the ocean did not respond locally to 
variability in this mode in spite of its being the most 
energetic ( 60% ) . Discussions in section 6 will consider 
some explanations for this result. 

A comparison between the hindcast ar) e and the ob- 
served 77 provides an assessment of the Markov model. 
Recall that ar} e is a time series at each grid point esti- 
mated solely from W E and two coefficients, ct and either 
* °[ ^ Jb e correlation between crq e and 77 is evaluated 
with CCA using three EOF modes from each dataset 
as input. (The EOF1 time series of each were actually 
quite similar, with a correlation of r = .73.) The re- 
spective CCA1 time series amplitudes (Fig. 14a) rep- 
resent the optimal correlation between the two fields. 

_ correspond quite closely with a correlation of r 
— .79. The respective CCA1 spatial amplitude contours 
are generally aligned and depict similar concentration 
?L\ anabl lty in the eastem portion of the gyre (Fig. 
cAci 7 CCA1 pattern is virtually identical to 77 
EOFl, and it is evident that the dominant mode of in- 
terdecadal dynamic height variability is largely ac- 
counted for with this simple model of direct Ekman 
pumping The CCA2 and CCA3 modes (not shown) 
had much smaller correlations (.41 and .22), and their 
spatial amplitude patterns did not portray the level of 
agreement as CCA 1 to suggest that local forcing of the 
same sort was explaining the remaining 77 variance. Of 
pa " note ’ a7 ]< did n °t produce a mode resembling 
77 EOF2, suggesting that an important fraction of the 
lowest frequency variability was not locally forced. 

The nature of the oceanic response to stochastic forc- 
mg is further illustrated in Fig. 15 with a simulation in 
which (4) is integrated over 100 years using computer- 
generated random white noise for W E , \ ~ l = l yr , and 
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a arbitrarily scaled. Decadal variability arises naturally 
from the random forcing. The lower half of the figure 
shows tj EOF1 to illustrate that the measured oceanic 
variability is similar in nature to the artificially forced 
autoregression model. Note that the time axes were de- 
liberately aligned to highlight the resemblance between 
the observed dynamic height and a similar segment of 
the random simulation. The similarity is purely coin- 
cidental, and this sequence was chosen for illustration 
from several such simulations. This provides a means 
to test the statistical significance of the CCA resu . 
Three random white noise sequences were pswsed 
through the autoregression model to generate artificial 
sequences like any 22.5-year segment of Fig. 15. The 


EOFs of these were supplied to a CCA with the first 
three tj EOFs, and the CCA1 coefficient was recorded. 
The process was repeated for 1000 trials to tabulate the 
CCA1 statistics. Using = 1 yrinthe autoregression, 
the largest coefficient was .78, or slightly less than .79 
obtained with arj e . Using X. -1 = 2 yr, 1.2% of the 
artificial coefficients exceeded .79. Given that \ 
generally fell between 0.5 and 2 yr with a median of 
1 yr, it is concluded that the CCA coefficient be- 
tween ct 7 ] e and tj is significant at about the 99% confi- 
dence level. 



6. Discussion 

The Markov model results remain deficient in some 
respects. As noted, the low-frequency variability as- 
sociated with tj EOF2 was not reproduced. This anom- 
aly is centered considerably westward of tj EOF1 Mid 
generally lags by a few years. Referring again to Fig. 
5, it is evident that the large-scale anomalies propa- 
gated generally westward, and it is common for two 
EOF modes to depict near quadrature phases of a prop- 
agating feature. In this case, it is a likely indication of 
long Rossby wave radiation from the eastern region oi 
local forcing. Interannual baroclinic Rossby waves 
have been documented along line P, near the soutiiern 
boundary of the present study area, by White and Ta- 
bata (1987). The conventional reduced gravity equa- 
tion for a forced long planetary wave in a L5 "^ a ^ r 
ocean, expressed in terms of dynamic height, is (Kelly 
et al. 1993) 


dt H dx g 


( 8 ) 


where (5 is the meridional planetary vorticity gradient, 
g ' reduced gravity, and local relative vorticity has been 





October 1995 


LAGERLOEF 


2253 



Fig. 1 1. Unfiltered W E EOF1 and W E EOF2 amplitude time series 
(light curves), which demonstrate the timescale separation between 
the climatic forcing term W E and the response term rj seen in Fig. 4a. 
The two-year running means (dark curves) are superimposed (see 

Fig. 10). 


ignored for large-scale L > L D . The internal Rossby 
radius L D is ~16 km based on the vertical eigenvalue 
solution to the mean buoyancy profile (see also Emery 
et al. 1984). The homogeneous solution to (8) is the 
nondispersive long Rossby wave with a phase speed C 
= Ld/3 105 km yr~ l at 55° latitude. The very long 
period 7 — 20 yr characterized by rj EOF1 and EOF2 
implies a wavelength L ~ 2000 km, or about 30° lon- 
gitude. This is the width of the study area and is con- 
sistent with the length scales of the features in these 
modes. Kelly et al. ( 1993 ) solved ( 8 ) with an idealized 
seasonally varying standing wave wind curl pattern to 
model wind-driven annual Rossby waves in the Cali- 
fornia Current. The solution contained standing and 
propagating modes with length scales dictated by the 
wavenumber of the atmospheric forcing and the free 
wave dispersion relation respectively. It is likely that 
similar dynamics are at work here on longer timescales, 
in particular the excitation of a long free wave with a 
wavelength determined by the timescale of the forcing 
and the dispersion relation rather than by the length 
scale of the forcing. The atmospheric wavelength as- 
sociated with the Aleutian Low exceeds 5000 km, so 
the system is not near resonance. A major discrepancy 
between the observations presented here and Kelly’s 
solution is that the standing wave response is out of 
phase with the forcing in the latter result. The low- 
frequency in-phase response shown here is explained 
by the damping term in (3). This suggests that a more 
general theory should include both damping and 
Rossby wave dynamics. Noting the similarity between 
( 1 ) and (8), this can be expressed as 

~ ^ = ~ aWE ~ Kr} ' (9) 



Fig. 12. Composite spectra formed by averaging the spectral den- 
sity of the lowest three EOF modes of rj (dark) and W E (light) re- 
spectively. Dashed curve is the idealized Markov spectrum (w 2 
+ A. 2 )"' with \ = 1 yr"' and uj = 2izf. 


where the damping embodies higher-order terms, 
which were ignored in the simplification of (8) from 
momentum and mass conservation equations. It is be- 
yond the scope of this work to attempt a rigorous der- 
ivation of such terms, but the empirical results provide 
numerical constraints, which may be useful in deter- 
mining appropriate scales. 

Another discrepancy with the Markov model results 
is the unbalanced oceanic response to the first two W E 
EOFs. As noted, W E EOF1 has a pattern centered over 
the gyre and has greater variance both in the higher and 



Fig. 1 3. Contours of the amplitude coefficient a derived from the 
Markov autoregression model. Note the eastward bias indicating 
dominant Ekman pumping response in the eastern region and weak 
response (a — 0) in the western portion. 








Latitude Amplitude (normalized) 
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otn, CCAl (dark) and i\ CCA1 (light) 



<xt\ c CCAl (dark) and t\ CCAl (Ught) 



Fig. 14. Canonical correlation analysis results between the ob- 
served dynamic height field tj and that estimated rj, from the observed 
Ekman pumping W E using the Markov stochastic forcing model to 
hindcast. (a) Temporal amplitudes 17 CCAl (light) and arj r CCAl 
(dark) with a canonical correlation of 0.79. The major events on 5- 
20-yr timescales are clearly replicated in both signals, (b) Spatial 
amplitude contours of 77 CCAl (light) and arj, CCAl (dark). The 
general agreement signifies that the CCAl time series share a similar 
spatial pattern of amplitudes between the observed and hindcast 
fields, particularly in the eastern portion. 


lower frequencies. Yet, there is little evidence of direct 
local response in accordance with the stochastic model. 
For some reason that is not clear, the ocean appears to 
have been more resonant with the spatial and temporal 
pattern associated with W E EOF2, which brought more 
vigorous Ekman pumping in the east and less in the 
center of the gyre. It is possible that the local midgyre 
response was obscured by the long propagating wave 
postulated above. An anonymous reviewer suggested 
that the observed response may be influenced by the 
eastern boundary current, bottom topography, or vari- 


ability associated with the Sitka eddy (Tabata 1982). 
Also, other forcing mechanisms such as surface heat 
flux may be important. The mean annual heat flux is 
near zero in the eastern part of the gyre and increases 
in magnitude toward the west according to new esti- 
mates by Tziperman and Bryan ( 1993 ) . A recent paper 
by Jacobs et al. ( 1994) describes a very long crested, 
long-lived Rossby wave propagating northwestward 
across the North Pacific subsequent to the 1982-83 El 
Nino. The northeastern extremity of the wave appears 
to have extended into the Gulf of Alaska according to 
GCM simulations, indicating a remotely forced oceanic 
teleconnection which may influence the dynamic to- 
pography on these timescales. 

To this point little has been said regarding the Sver- 
drup transport. A barotropic Sverdrup balance has been 
described in deep North Pacific currents on shorter than 
yearly timescales (Koblinsky et al. 1989; Niiler and 
Koblinsky 1985). Regarding the baroclinic mode, it is 
commonly accepted that a quasi-steady Sverdrup bal- 
ance will account for the interannual large-scale vari- 
ations in the meridional transport and the Alaskan 
Stream. It was noted above, however, that prior obser- 
vations do not support this on seasonal to interannual 
timescales. Wunsch and Roemmich (1985) also ques- 
tioned the basic premise in the subtropical North At- 
lantic. In the quasi-steady version of ( 8 ), a negative 
W E anomaly pattern over the eastern portion should 
correlate with decreased drj/dx. This implies an 77 
anomaly with a relative high at the gyre center and low 



Fig. 15. A simulation of stochastically forced variability using the 
Markov model (Eq. (4)] with computer-generated white noise for W E 
and \ = 1 yr _t (top panel). The low-frequency variations produced 
are similar in nature to observed 77 variability, as demonstrated with 
j; EOF1 (bottom panel). Note that the time axes were deliberately 
aligned to highlight the similarity between the observed dynamic 
height and a segment of the random simulation. The similarity is 
purely coincidental, and this sequence was chosen for illustration 
from several such simulations. 
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at the eastern boundary rather than the observed pat- 
tern. To consider whether the quasi-steady assump- 
tion is valid, note that the terms on the left side of 
(8) are similar in magnitude for decadal timescales of 
T - Lp~'Lo 2 « 10-20 yr( with L « 1000-2000 km). 
This is the period of the long Rossby wave that would 
bring about the geostrophic adjustment to a large-scale 
change in the wind curl. Quasi-steady Sverdrup balance 
would require the wind timescale to be many decades 
before the first term in Eq. ( 8 ) becomes negligible. This 
slow adjustment was also described by Liu ( 1993 ) with 
an ideal two-layer subtropical gyre model response to 
a single step change in Ekman pumping. The timescale 
required for Sverdrup adjustment remains much longer 
than the duration of this 22.5-year record. Conse- 
quently, it should be no surprise that the baroclinic re- 
sponse seen here is not consistent with theoretical Sver- 
drup transport. The transient response is a local Ekman 
pumping effect to which the meridional flow adjusts 
geostrophically rather than responding directly to the 
wind curl in Sverdrup balance. 

Randomly forced decadal variability as seen in Fig. 
15 will obscure the signature of weak 18.6 year nodal 
tide SST variation posed by Royer (1993), necessi- 
tating very long and stable records to verify the har- 
monic constituent. Royer reviewed several examples 
of — 18-22 year variability in climate indices in North 
America, the North Atlantic, in various fish stocks, 
and other cases. A significant spectral peak in this 
band for the Sitka, Alaska, coastal air temperatures, 
which are coherent with adjacent SST, strengthened 
the argument in favor of the nodal tide. However, it 
remains unresolved whether this variability originates 
in the ocean (i.e., the nodal tide) or in the atmosphere. 
The results here demonstrate that low-frequency pe- 
riodicity in the atmospheric circulation would be 
transferred to the ocean and suggests that the stochas- 
tic forcing alone accounts for the observed ocean vari- 
ability. 

It has been noted that the climate changes in the 
North Pacific during this period included deepening of 
the Aleutian low, southward displacement of Pacific 
storm tracks, lower SST, and a deeper mixed layer 
across the central North Pacific and higher SST along 
the North American coast and Alaska. These are all 
consistent with the observations presented here. The 
cause of this large-scale interdecadal anomaly appears 
to be related to tropical Pacific SST (Graham et al. 
1994). However, the point to be made here is that cli- 
matic variability witnessed in the Alaska Gyre is clearly 
a manifestation of much larger scale happenings in the 
North Pacific ocean -atmosphere climate interactions. 
Trenberth and Hurrell ( 1994) presented the large-scale 
average NP pressure index as a robust gauge of this 
climatic signal. A smoothed contour plot by year and 
month (their Fig. 2) shows an interdecadal oscillation 
of midwinter SLP with minima in the mid- 1920s, early 
1940s, early 1960s, and mid-1980s. These coincide al- 


most precisely with the peak Alaska coastal air tem- 
peratures shown by Royer ( 1989, his Fig. 8), suggest- 
ing these air temperatures are the result of warm air 
advection by southerly wind anomalies associated with 
the Aleutian low. This bears on the discussion of the 
nodal tide above and lends added support to the hy- 
pothesis that the interdecadal variations in the Gulf of 
Alaska are related to the large-scale variations in at- 
mospheric circulation. 

7, Summary and conclusions 

Circulation in the Alaska gyre is shown to have 
large-scale variations over decadal timescales during 
the period 1968-1990. These are related to previously 
documented interdecadal SST variations, with wanner 
SST associated with increasing dynamic height in the 
north and eastern parts of the gyre. The gyre center was 
farther east and circulation appeared stronger during 
the cooler SST decade of the 1970s than the warmer 
SST decade of the 1980s. The Aleutian low (PNA in- 
dex) intensified during the interim, and the associated 
wind stress anomalies forced a slowly varying dynamic 
height anomaly across the northern and eastern part of 
the gyre through Ekman convergence. This had the ef- 
fect of displacing the gyre’s low somewhat to the 
WSW, with an apparent relaxation of circulation in the 
central and eastern gyre in the 1980s rather than a gyre 
spinup as might be anticipated with the intensified 
Aleutian low. Reed (1984) had suggested that accu- 
mulated pycnocline displacement associated with local 
wind stress curl accounted for the apparent westward 
displacement of baroclinic flow observed between two 
surveys in the early 1980s, and the present results sup- 
port that hypothesis. The relationship between patterns 
of wind stress and the dynamic height were not con- 
sistent with poleward transport response according to 
the Sverdrup balance. The data instead support dy- 
namic height response through vertical pycnocline 
forcing by Ekman pumping and circulation changes by 
geostrophic adjustment. 

The wind curl spectrum was white at periods > 100 
days. The transient baroclinic ocean response had a 
damping timescale of ~1 yr. By approximating the 
forced climatic response in terms of a first-order Mar- 
kov autoregression process, the development of deca- 
dal-scale variations was explained as the low-frequency 
oceanic response to stochastic wind forcing. Using this 
model with observed winds, hindcast dynamic height 
temporal and spatial variations were consistent with 
those observed. A simulation using artificial white 
noise forcing shows the evolution decadal variations 
similar in nature to the observed variations. These re- 
sults demonstrate that much of the observed dynamic 
height and SST interdecadal variability can be attrib- 
uted to stochastic wind forcing. 
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APPENDIX 

Data Distribution and Objective Analysis 

The geographic distribution of observations used in 
the dynamic height analysis is shown in Figs. A1 and 
A2. These consist of 9382 XBTs and 2769 STD/Nan- 
sen casts, for a total of 12 151 profiles in the study area 
spanning the time period from 1968 to mid- 1990. The 
latter consist of traditional hydrocasts using sample 
bottles and reversing thermometers as well as the mod- 
em CTD profiles. The method for estimating dynamic 
height from XBT data is given by Lagerloef (1994). 
The point observations of 0-450 dbar dynamic height. 
rj were interpolated to a fixed grid at the intersection of 
each whole degree of latitude and longitude. Note that 
this scale does not resolve the details of the shelf and 
coast; thus contours may appear at times to cross the 
shelf and the superimposed shoreline as an artifact of 
the coarse grid. Only oceanographic station depths ex- 
tending to >450 m were used in the analysis, so data 
remain representative of oceanic and not shelf mea- 
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Fig. Al. Geographic distribution of the 9382 XBT stations, 
1968 to mid- 1990. 


2769 STD/Nansen Stations: 1968 through mid- 1990 



Fig. A2. Geographic distribution of the 2769 STD/Nansen stations, 
1968 to mid- 1990. 


surements. The temporal resolution of the interpolated 
fields was set at 16 samples per year so that the Nyquist 
frequency would be sufficiently high to resolve annual 
and semiannual signals should they be of interest. 

The interpolation followed common practices in sta- 
tistical objective analysis (e.g., Thiebaux and Pedder 
1987). This required the specification of a correlation 
model, and a simple Gaussian form was used: 




. (Al) 


where A*, Ay, and Af are separations in the east, 
north, and time dimensions, respectively; and S x , S y , S, 
are the respective decorrelation scales. The mean field 
was estimated by first taking the ensemble mean of all 
observations as a first guess and subtracting it from 
each observation. The interpolation of the residuals was 
made to each grid point using only the first two terms 
of (Al) with large scales S z - 30° longitude and S y 
= 10° latitude to ensure a smooth field. Adding back 
the ensemble mean gave the gridded mean field shown 
in Fig. 1 . The space-time objective analysis began with 
an interpolation of the mean field to each observing 
location and subtraction to yield the residual from the 
local mean at that location. The residuals were then 
interpolated to the space-time grid using scales S x = S y 
= 200 km and S, = 0.25 yr to produce the dataset for 
subsequent analysis. In order to avoid the inversion of 
very large correlation matrices, an influence radius was 
employed to screen observations outside the neighbor- 
hood of the interpolation grid point. For computing the 
mean field,, this was arbitrarily set to r — [r; + r y ] 1/2 
with r x — r v = 100 km. For the subsequent space-time 
interpolation, r — [r? + rj 4- rj] 1/2 with r x = r y — 300 
km and r, — 0.375 yr. These radius components were 




■MjagaaasgaaEMgfflaE^ 


'i-i-V-ai-.. . .- .-: -1.. ...v 


October 1995 LAGERLOEF 2257 


1.5 times their respective decorrelation scales, so that 
observations with a correlations less than about 
exp[— (1.5) 2 ] ~ 0.1 were ignored. In general, the 
number of observations included in each interpolation 
point ranged between about 50 and 200. Interpolation 
of the 2 deg X 2 deg monthly COADS wind stress was 
computed with scales S x — S y = 1000 km and S, = 0.25 
years to filter the higher wavenumber errors and reduce 
noise in the subsequent finite-difference computation 
of the wind curl. TTie coarse 2 deg X 2 deg COADS 
grid implies that the coastal winds and the related 
coastal Ekman convergence are not well resolved. 
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